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Abstract 

Within framework of the quantum calculus, we represent the partition function and 
the mass exponent of a multifractal, as well as the average of random variables 
distributed over self-similar set, on the basis of the deformed expansion in powers 
of the difference q — 1. For the partition function, such expansion is shown to be 
determined by binomial-type combinations of the Tsallis entropies related to man- 
ifold deformations, while the mass exponent expansion generalizes known relation 
t q = D q (q — 1). We find the physical average related to the escort probability in 
terms of the deformed expansion as well. It is demonstrated the mass exponent can 
acquire a singularity that relates to a phase transition of the multifractal set in the 
course of its deformation. 
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1 Introduction 



Fractal conception [1] has become widespread idea of contemporary science 
(see Refs. [2] - [5], for review). Characteristic feature of fractal sets is known 
to be the self-similarity: if one takes a part of the whole set, it looks like the 
original set after a deformation A ^ 1. Formal basis of the self-similarity is the 
power-law function F ~ X h with a Hurst exponent h (for time series, value 
F is reduced to the fluctuation amplitude and A is the interval size within 
which this amplitude is determined). While the simple case of monofractal is 
characterized by a single exponent h, a multifractal system is described by 
a continuous spectrum of exponents (singularity spectrum) h(q) [6]. Along 
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this way, a self-similarity degree q represents the exponent of a homogeneous 
function being characteristic function of self-similar systems [7] (so, within 
nonextensive thermostatistics, this exponent expresses the escort probability 
Pi oc p\ throughout the original one Pi [8,9]). In physical applications, a central 
role is played by the partition function Z q oc l T ^ (I is a characteristic size 
of boxes covering multifractal) to be determined by the mass exponent r(q) 
connected with the generalized Hurst exponent h(q) by the relation r(q) = 
qh(q) - 1. 

Since fractals are scale invariant sets, it is natural to apply the quantum 
calculus to description of multifractals because quantum analysis is based 
on the Jackson derivative which yields variation of a function with respect 
to the deformation of its argument [10,11]. This Letter is devoted to such a 
description with using deformed series of expansions of the partition function, 
the mass exponent and the averages of random variables over deformed powers 
of the difference q — 1. As shows the consideration in Section 2, the coefficients 
of the partition function expansion are reduced to binomial-type combinations 
of the Tsallis entropies related to manifold deformations A m with growing 
powers m = 0, 1, . . . Respectively, Section 3 shows that escort probability is 
generated by the mass exponent whose expansion into deformed power series 
represents the generalization of well-known relation r q = D q (q — 1) with fractal 
dimension spectra D q corresponding to manifold deformations. According to 
Section 4, above expansion permits also to expresses physical averages over 
the escort probability Pi in terms of related averages based on the ordinary 
one Pi. Section 5 contains discussion of results obtained, and in Appendix we 
adduce necessary information from the quantum calculus. 



2 Partition function 

Following the standard scheme [6,2], we consider a multifractal set covered 
by boxes % = 1,2, ... ,W with number W — > oo. The peculiarity of self- 
similar sets is that the probability to occupy a box % is determined by its size 
li — > according to the power-law relation pi = If characterized by a Holder 
exponent a > 0. All properties of a multifractal is known to be determined by 
the partition function 

w 

Z q = Z q {p t }:=Y / Pl (1) 
i=i 

which according to the normalization condition takes the value Z q — 1 at 
q — 1. Because pi < 1 for all boxes i, the function (1) decreases monotonically 
from maximum magnitude Z q = W related to q — 0. In the simplest case of 
the flat distribution pi = 1/W fixed by the statistical weight W ^> 1, one has 
the exponential decay Z q = W l ~ q . 
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Our approach is based on the expansion of the partition function (1) into the 
deformed series (23) in powers of the difference q — 1 [10] 

oo q( n ) 

n=0 L' 4 M- 

Here, the deformed powers (q— 1)^ are determined by the binomial (24), while 
the kernels = S^{pi} are defined by the n-fold action of the Jackson 
derivative: 

S^ = -(qV^ n Z q \^. (3) 

By using the method of induction with accounting the definitions (22) and 
(1), one obtains the explicit expression 

= -(\-l)^£(-l)™(")z xm (4) 

m=0 \ m / 

where ( n ) = —rr^ — tt are binomial coefficients. In absence of deformation 

\mj m\{n—m)\ 

(A = 1), all coefficients (4) equal zero, apart from the single term S[°^ — — 1. 
Inserting these coefficients into the series (2) gives Z\ — 1 as desired. 

It is easily to see that the set of kernels (4) with n > is expressed in terms 
of the Tsallis entropy [8] 

c r -i YsiPi — 1 Z\ — l 

With growth of the deformation parameter A, this entropy decreases mono- 
tonically with inflection in the point A = 1 where the expression (5) takes 
the Boltzmann-Gibbs form Si = — X^Pi m (Pj)- m the case n — 0, the coeffi- 
cient (4) does not depend on the deformation parameter A to give the value 
S| 0) = -1. At n = 1, one obtains [12] 



S?{ Pl }:=-^4 1 ="^ = W (6) 



A- 1 



q=l 



A - 1 



Respectively, in the second order n — 2, Eq.(4) yields 



With accounting the definition (5) one has 
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In arbitrary order n, combination of equations (4) and (5) arrives at the sum 

W (A — 1)™ 

which expresses in explicit form a generalized entropy of a multifractal by 
means of contributions of the Tsallis entropies related to different powers of 
the set deformation. In the limit A — > 0, when S\ ~ W, accounting of all terms 
in the sum over m = 0,1, ... ,n arrives at cancellation of the denominator 
(A - 1)™ so that S<5 n) -> W; similarly, at A -> 1, one obtains the Boltzmann- 
Gibbs limit — > Si = — X^Pihife) as desired. Finally, in the case A — > oo, 
where S\m ~ A~ m , the term related to m = is reduced to zero so that 
main contribution gives the term with m — 1 to arrive at the sign-changing 
asymptotics ~ (-l) n+1 A~ n . 

For the flat distribution, when Z\ = W l ~ x , the dependence (5) is characterized 
by the asymptotics S\ ~ W [1 — ln(W)A] in the limit A 1 and S\ ~ 1/A at 
A > 1. As a result, with the A growth the coefficients (4) increase from the 
value Sq^ = W to the Boltzmann magnitude S[ n ^ = W ln(W) and then decay 
to the asymptotics ~ (— l) n+1 A~ n . 



3 Characteristic function of escort probability 



Following pseudo-thermodynamic picture of multifractal sets [9] , let us define 
effective values of the free energy r q , the internal energy a and the entropy /: 

_ HZ q ) EjPjhm f _ E l P i lnP l 

Tq - ln(Z) ' a ln(l/0 ' 1 - ln(l//) " 1 UJ 

Here, / 1 is the size of characteristic box in phase space, Pi and Pj are 
original and escort probabilities to be connected with the definition 

p ■= — q = #• (11) 



E t pJ z, 



Inserting this equation into the second expression (10), one obtains the well- 
known Legendre transformation [6,2] 

r q = qa q - f(a q ) (12) 

where q plays the role of the inverse temperature and the internal energy is 
specified with the state equation 

a a = (!3) 
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It is easily to convince the escort probability (11) is generated by the charac- 
teristic function being the mass exponent 



_ ln(Zg) _ InCEiP?) 
q -~ ln(Z) ~ ln(/) ' 1 } 



Really, one has 



P- - ln(/) v dT « - g^ dln ^ (15) 
q dpi om(pi) 

By analogy with Eq.(2), one can expand the function (14) into the deformed 
series 



x n ( n ) 



^SH^- 1)r (16) 

being the generalization of the known relation r q = D q (q — 1) connecting the 
mass exponent r q with the multifractal dimension spectrum D q [2]. Similarly 
to the equations (3) and (4), the kernel is determined by the n-fold action 
of the Jackson derivative: 



4 n) := (<^HA-irEHr(>». (17) 



m=0 



Noteworthy, the term with n = is absent in the series (16). At n = 1, one 
reproduces the ordinary relation = t\/(A — 1), while the kernels 
with n > 2 include terms being proportional to T\m/ (A — l) n to correspond to 
manifold deformations A m , 1 < m < n. Generally, the definition (17) yields a 
hierarchy of the multifractal dimension spectra related to different multiplica- 
tion factors n of the set deformation. 

In trivial case of the flat distributional], when the partition function is Z\ = 
W 1 ~ x , the definition (14) yields the mass exponent T\ = D(l — A) with the 



fractal dimension D = which tends to D — 1 when the size of covering 

boxes / goes to the inverse statistical weight 1/W. Thus, we can conclude the 
flat distribution relates to a monofractal with dimension D\ = y~^y tending 
to a smooth one-dimensional set at deformation A — * 1/7. 



1 It is worthwhile to stress the expression (15) is not applicable in this case because 
fixation of the statistical weight W makes impossible variation of the probability 

Pi = i/w. 
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4 Random variable distributed over multifractal 

Let us consider an observable fa distributed over a multifracal set. According 
to Ref. [8], its mean value 

(fa:=J2faPi (18) 

i 

is determined by the escort probability (11). With accounting Eqs. (15) - (17), 
the average (18) can be expressed in terms of the sum 

where the specific average 

(0> A ™ = ]>>P; (A™) (20) 

i 

is introduced. At m — 0, the escort distribution Pj (A m ) is reduced to the 
original one p; so that related terms in the sum (19) are proportional to the 
ordinary average (fa) l = J2i <PiPi- On the other hand, in absence of deformation 
(A = 1), the only term with n = determines the mean value (18): {fa) = (fa) v 
In the limit A — * 0, the escort probability reduced to flat distribution Pj(0) = 
1/W for arbitrary form of the ordinary distribution so that average (18) is 
determined by the simple expression (fa = W^Y.ifa- Finally, in the case 
A — > oo, main contribution is given by the terms with m = and m = 1 so 
that the average (19) takes the form 

(<f>) — (0)1 + (q- i)((0)a _ (21) 



5 Conclusions 

A principle peculiarity of above consideration is that the expansions (2), (16) 
and (19) depend of both deformation parameter A and self-similarity degree 
q. As shows the example of using the flat distribution in the end of Section 3, 
the former determines the fractal dimension D\ = j^-j^y in the simple case of 
the monofractal generated by deformation A ^> 1 of initial cell of size I <^ 1. 
With passage to general case of multifractal, we obtain the fractal dimension 
spectrum defined by the series (17). Along this line, the self-similarity degree 
q plays the role of a free parameter whose variation describes the multifractal 
spectrum, while value of the deformation parameter A fixes the scale of the 
multifractal resolution to be determined by external conditions in the case of 
natural objects. 
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As pointed out in Introduction, our study is based on the Jackson deriva- 
tive (22) which gives variation of a function with respect to scale choice of 
its argument. This derivative expresses the coefficients of the deformed Tay- 
lor expansion (23) in powers of the deformed difference q — 1 defined by the 
binomial (24). In this way, the expansion (2) of the partition function (1) 
is determined by the coefficients (9) being binomial-type combinations of the 
Tsallis entropies (5) related to manifold deformations A m with growing powers 
m — 0, 1, . . . According to definitions (10) and (11) the mass exponent (14) is 
expressed by the Legendre transformation (12) to generate the escort probabil- 
ity (15). On the other hand, the expansion of the mass exponent into deformed 
power series (16) represents generalization of known relation r q = D q (q — 1) 
with fractal dimension spectra (17) corresponding to manifold deformations. 
Moreover, making use of the deformed expansion (23) arrives at the expression 
(19) of the physical averages (18) corresponding to the escort probability P { 
in terms of the averages (20) related to multiple deformations. 

Under deformation A, above multifractal characteristics vary in the following 
manner. The generalized entropies (9) tend to the statistical weight W in the 
limit A — > 0, take the Boltzmann-Gibbs magnitude S[ n ^ = —J2iPi^ n {Pi) at 
A = 1 and arrive at the sign-changing asymptotics ~ (— l) n+1 A~ n when 
A — > oo. What about the relation (19) between the physical average (<fi) = 
^Zi<piPi and the ordinary one = J2i<friPi, one has the simple expression 
(<f)) = W~ x J2i<fri i n the limit A — > 0, the trivial relation (<fi) = (<fi) 1 in absence 
of deformation (A = 1), and the difference ((f)) — ((f)) 1 being proportional to 
the factor q — 1 at A — > oo. 

From physical point of view, the most tractable is the expansion of the mass 
exponent (16) whose coefficients represent fractal dimension spectra. As shows 
the equality (17), the fractal dimensions related to different powers n 
of the series (16) can change their signs with variation of the deformation 
A. Because, within the pseudo-thermodynamic formalism, the mass exponent 
r and the self-similarity degree q relate to the free energy and the inverse 
temperature respectively, such a change means that thermodynamic potential 
can obtain a singularity at some deformations A. Within framework of the 
physical presentation [9], this singularity means the phase transition in the 
course of the deformation of the multifractal set. 



Appendix 



Quantum analysis is known to be based on the Jackson derivative [10,11] 
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where A is a deformation parameter. The deformed Taylor expansion reads: 

(x - a)S 2) 



/(*) = /(a) + ^TT# ^/tt 



[1]; 



+ 



[2]a! 



IX. 



A (^)7( 

-s^ (*)•«■> 

Here, the deformed binomial 

(x + yt ] = (x + y)(x + Xy)(x + \ 2 y) ...(x + X^y) 



+ 



(23) 



= E 

m=0 



m 



m(m— 1) 



(24) 



is determined by the coefficients 



m 



[m\ x \[n - m] A ! 



(25) 



where generalized factorials [n]\\ = [1]a[2]a • • • Na are given by the basic de- 
formed numbers 

\n 1 

Ma = p ( 26 ) 
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